NLO corrections to timelike, spacelike and double deeply virtual Compton scattering. 
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We calculate the 0[as) corrections to the timelike, spacelike and double deeply virtual Compton 
scattering amplitudes in the generalized Bjorken scaling region. Special attention is devoted to 
studies of the difTerence between the next to leading order timelike and spacelike coefficient functions, 
which plays for this process a role analogous to the large K factor which was much discussed in the 
analysis of inclusive Drell Yan cross sections. Also in the present studies the timelike nature of the 
hard scale gives rise to new absorptive part of the amplitude and to the presence of characteristic 
TT^ terms, which can potentially lead to sizeable corrections. 

PACS numbers: 13.60.Fz , 13.90.+i 



INTRODUCTION. 



I Data on deeply virtual Compton scattering (DVCS) are now available from various experimental settings [l| and 
different strategies are proposed to extract from them the physical knowledge on nucleon structure encoded in 
generalized parton distributions (GPD) 0, 9- These attempts are usually based on a leading order QCD analysis, 
although the importance of next order terms has often been emphasized, in particular with respect to the dangerous 



^ '' factorization scale choice dependence @ . Historically, one can note that the understanding of inclusive reactions (Drell 
Yan reactions, large pr hadron or jet production) in the framework of collinear QCD factorization has waited for an 
^ • analysis including next-to-leading order (NLO) (or even next-to-next-to-leading order) corrections. Indeed, complete 
I I \ NLO calculations are available for the DVCS reaction and there is no indication that they are negligible in the 
kinematics relevant for current or near future experiments. Deeply virtual Compton scattering is only one case of the 
general double DVCS(DDVCS) reaction 

l*{q„r)N ^^*{qout)N' , (1) 

where the final photon is on shell, g^^j = 0. The converse case where g|„ =0, often called timelike Compton scattering 
(TCS), has been theoretically discussed at medium [l^ and very large [ll| energy and first data are being analyzed 
hj. The double DVCS case has been discussed in Ref. [l3| . 

It has been shown that the understanding of DVCS data needs higher order calculations for a reasonable extraction 
of GPDs to be possible This is likely to be even more the case for TCS. Indeed, TCS and DVCS amplitudes are 
identical (up to a complex conjugation) at lowest order in as but differ at next to leading order, in particular because 

_ ^ ^ of the quite different analytic structure of these reactions. Indeed the production of a timelikc photon enables the 

, production of intermediate states in some channels which were kinematically forbidden in he DVCS case. This opens 
\^ • the way to new absorptive parts of the amplitude. Soon, experiments will be performed at JLab at 12 GeV which will 
5^ I enable to test the universality of GPDs extracted from DVCS and from TCS, provided NLO corrections are taken 
into account. Experiments at higher energies, e.g. in ultraperipheral collisions at RHIC and LHC, may even become 
sensitive to gluon GPDs which enter the amplitude only at NLO level. 

Former experience with inclusive deep reactions teaches us that NLO corrections are likely to be more important in 
timelike reactions than in corresponding spacelike ones. The well-known example of the Drell Yan K-factor teaches us 
that NLO corrections are sizeable in timelike processes, because of in factors coming from log{~Q^ / fijp) terms which 
often exponentiate when soft gluon resummation is taken care of 0,[ll. 

The results for TCS should be indicative of other exclusive reactions with a timelike scale, as TriV l^l~N' 
discussed in [l7| which may be accessed in the Compass experiment at CERN or at J-Parc, e^e~ — > ^tttt discussed 
in fisl to be compared to 7*7 — !■ tttt analyzed in [l^, or 7*A^ — > A^'Trfg^] to be compared to NN' 7*7r analyzed in 

Our calculations are performed along the lines of Ref.0 (see also Q)- Those earlier results where obtained in an 
unphysical region of parameters space and then by analytical continuation (due to simple analytical structure of hard 
DVCS amplitude) extrapolated to the physical region of DVCS. Not restricting the parameters enables us to get the 
full result for the general kinematics (including TCS, DVCS and DDVCS). In earlier analysis the factorization scale 
fip was kept equal to the hard scale Q^. In our calculation they are independent, so one can check factorization scale 
dependence. We calculate only the symmetric part of the amplitude which is dominant for the phenomenological 



2 



analysis, as the main features of scattering amplitudes of DDVCS, DVCS and TCS are already clearly seen. We 
simplify the kinematics by restricting ourselves to the forward (t — imm) region. We leave the phenomenological 
analysis of our results to a future publication. 



II. PRELIMINARIES. 




Figure 1: Incoming and outgoing photon virtualities as a function of a;_B. To get the physically interesting case in which 
incoming photon (solid line) is spacelike and outgoing photon (dashed line) is timelike, one has to choose > for xb > 
and < for xb < 0. DVCS kinematics corresponds to a;s = ^, and TCS to xb = — 



As in [6| we describe the kinematics of general Compton scattering in a symmetric way i.e. in the Bjorken limit, 
and in the forward limit {q^^t — 0, ~ 0) momenta are assigned as follows: incoming photon qin = (q — ^p), outgoing 
photon {qout = 9 + S,p), incoming proton P = (1 + and outgoing proton P' = (1 — ^)p, where: 



P 



q = 



so: pn = 1, s = [p+qY = 



-Q^ and xb 



s+C 



0,0,1), 
^(1,0,0,-1), 

-XBP+ n, 

2xB 

We use the following vector decomposition k" = k'^p" 



(2) 



' + k° 



With above definitions we arrive at the following equations for incoming and outgoing photon virtualities: 



2 



Xb 



2 

qout 



Xb 



(3) 



From this general kinematics we can get as a limit some physically interesting cases. It is easy to check that to get 
incoming photon momentum spacelike and outgoing photon timelike, one has to choose > for > and < 
for Xb < 0. Deeply virtual Compton scattering is restored for xb = C ^^'^ > 0; timclike Compton scattering 



for a;_B = and Q 



-Q' < 0, and double deeply virtual Compton scattering for < < ^ and Q > or 



> Xb > — C and < 0. This is illustrated by Fig. [T] which shows incoming and outgoing photon virtualities as a 
function of xb. 

We perform our calculation in MS scheme, with D = 4 + e regularizing infrared divergences, as all ultraviolet 
divergences cancel out. In the following we shall study only the symmetric part of the full Compton scattering 
amplitude since it is phenomenologically the dominant part. Its factorized form reads: 



dx 



J2 f'>{x)F'>{x) + f\x)Fa{x) 



(4) 
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Renormalizcd GPD's arc defined as in [6|, by: 



1 r dx 
— f 

27r 



P' 



P )n 



(5) 



The connection between the bare quantities Fq, Fg and the renormalizcd ones in MS is given by: 



Fi{x) = F-i{x) 



F3{x) = FS{x) 



1 1 . einl 
e 2"' 4^ J 

- + J In K^'ix, x') ® F^ix') , 



ln-i4^ «)P9(x') 



(6) 



where evolution kernels K"^"^, K'^^, K^^ , K^'' may be read from [j], (g) stands for integration over common variable from 
-1 to 1. At the NLO of the process studied in this paper the parts with K^^ and K^^ do not contribute, since the 
gluon contribution is of the 0{as)- 

In Eq. ^ unrenormalized coefficient functions contain infrared divergencies, and are given by: 



CI). 



47r^^ 



^coU + "-'1 



(7) 



is calculated using the following relation with 57 — s- (77 hard scattering amplitude M^'^, given by diagrams shown 
on Figs. [21 [3] and m without attachement of external spinors of the t-channel quarks 



2"? 



2^1— Tr 

D-2 



(8) 



In dH]) the factor -^^^ is a part of the projection operator in Lorentz indices on the symmetric part of the full Compton 
scattering amplitude in Eq. ([4]). Factor 2 is related to the definition of quark F"^ given by formula ([5]). 

is calculated using the following relation with 37 — s- 97 hard scattering amplitude j\4'^'^°'^ , given by diagrams 
shown on Figs. and [H] without attachement of external polarization vectors of the t-channel gluons 



1 



-2x 



2 {x — xb + ie){x + xb — ie) {D — 2) 



{D~2) 



(9) 



Similarly to the quark case, factors j^z2 ^^d (j)L2) ^''^ parts of the projector operators on symmetric two photon and 
two gluon states, respectively. The factor ^ is the combinatorial factor which appears due to the condition that on the 
gluonic target we reproduce usual contribution of six diagrams shown in Figs. [5]and[6l The factor g-\^i^){^x+x b -ie) 
requires more explanations. It appears since in the axial gauge n ■ A = Q we have the relation: 



P' 



-2x 



A 



P )5T„ 



(x - XB + ie){x + XB - ie) 



P' 



Si 



n F, 



n 

2 



/3 • 



(10) 



The structure of denominators in (|10p is not fixed by the gauge condition n ■ A = alone. This arbitrariness is 
due to the presence of the residual gauge. It is fixed by additional boundary conditions involved in the factorization 
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procedure of the whole scattering amplitude of the given process. Here we fix it in agreement with the structure of 
denominators in the quark Born coefficient function for general double DVCS kinematics: 



^0 (DDVCS) ^ 



1 



1 



(11) 

— xb + i£ X + xb ^ 

In particular, in the case of DVCS where xb = ^ we obtain the standard expression (a; — ^ + ie)(x + ^ — ie) (see Q)- 
In the case of the TCS where xb = —S, this product becomes (.t + ^ + ie){x — ^ — ie). Detailed calculation of T'' and 
will be presented in the section UIIl 
If the following relations between Born coefficient function, infrared divergent terms and evolution kernels hold: 

Clui^') = Cl{x)®K'i^{xy). 
one can rewrite the full amplitude in the fully factorized form: 



(12) 



A^"" = 5^'' j dx 



n F 

E 

. 9 



Ti{x)F'i{x)+TS{x)F3{x) 



(13) 



where renormalized coefficient functions are given by: 



^coll ' 



T^^Cf + lln(^ 



(14) 



In the next section we will describe one-loop calculations necessary to obtain the above coefficient functions, in more 
details, as they can be useful in the calculations of similar processes (for example [l7|). 

III. INTEGRALS 

A. Integrals with two propagators. 

We start with a detailed description of the diagram shown on Fig. [2l Although this calculation is very simple, it 
reveals some characteristic features of the full calculation, and some pattern of the analytical structure of the result. 






Figure 2: Self energy correction to 57 — !> 57 scattering amplitude 



The symmetric part of the amplitude is given by: 

1 



Tr [M^" 



(dk)- 

[{q + x-pf + ie]2 }^ ' {(k + q + xpY + ie] [fc^ + ie] 



(15) 



where (dfc) = /i 



4-_D d"k 



ir, and Cp 



We have two types of integrals to perform: 



60 



{dk)- 



[{k + q + xpY + ie][/c2 + ie] 
^'^^\{k + q + xJ^^\^eW+^e] = 'V" + "^^^^^ ' 



(16) 
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bo may be shown to be equal (pay attention to the difference between e and e): 

1 



. ^_ - j.( ^A r(i + f)r(i + f) f XB- 

" (4^)2 (47r^2)f V 2/ r(2 + e) \^ xb 



(17) 



When we add the diagram with external photon lines crossed, which is given by the xb ^ ^xb substitution, we get 
the following result for the sum of those two diagrams with the self-energy corrections: 



-9t 



1 



47r {Atth 



2\^ 



1 



XB [Q^^+ie] 
1 

XB [-Q'^+^s] 



Q' 



,Xb -X 



— le 



which in the DVCS and TCS limits results in: 



Tr [M'^" 



'DVCS 



-9t 



e^asCp ( 



in 



'9t 



477 



[x 




ie] 




1 




[x 


-e + 


ie] 




1 




[x 




ie] 




1 




[x 




ie] 



2 

e 

2 

e 

2 
e 

2 
e 



Xb 

^2 Xb + X 

Xb 



— le 



(18) 



1 - log \\ + 
1 - log (^1 - 
l-log(^-l 
1 - log {-\ 



X 






iej 


X 






ie^ 


a; 






— ie 


X 






— ie 



(19) 

We notice that in the TCS case we have £, + ie contrarily to the £^ — ie present in the DVCS. There is also an overall 
minus sign under the logarithm, coming from the different sign of Q^. 

B. Integrals with three propagators. 

In this section we will describe in a detailed way the calculation of the diagram shown on Fig. [3l because all other 
diagrams with three, and some of the diagrams with four propagators in the loop, may be calculated in a similar way. 
The symmetric part of the amplitude with the right vertex correction is given by: 




Figure 3: Right vertex correction to qj — >■ qy scattering amplitude 



(q + xpY + ie J [{k + g)2 + ie] [{k — ^p)^ + ie] [(fc — xp)'^ + ie] 

(20) 



We start with the integration over k . There are three poles placed at: 

fc2 — ie 



k; 



2{y - x) 



6 



k2 - 



le 



l-^] -ie 

Xb ' 



2(2/ -xs) 



(21) 



For various values of y we close the contours of integration in the upper or lower half plane, in such a way that we 
avoid catching the pole. Irrespectively of the ordering of .t, xb and ^, we arrive at: 



1 



(q + x'pY' + 



; £ % I (dk^) Res,-/ - ^P^l Res,-/ 



where {dkj_) = /i ^-"^"^ V ^ k^ ■ Residua of the first and the second pole are given by: 



Res,-/ = 
Res,-/ = 



?/ - a; 



2(x-XB)(a;-0 [A;2^] 



kl^Q 



2 (£b_leHi^:£) 

"J- ^ xb(x-xb) 
"2 + /32fci 



jL2 _ ^2 (^B-y)fa-g) _ „> 



2(C-a;s)(C-2:) [^2] 
where a, and /?i are defined by the value of the numerator at the correspondent pole 

After wc perform the integration over k~ we arrive to: 

1 



{q + xpY + i£ 



? / idk.) 



27r 



y~x 



ai + f5ik]_ 



2{x-xb){x-^) ^j^2j j„2^ Q2 {^B-y)(y-x) 



dy 
2tt 



(dk^) 



y-i 



xb(x~xb) 
"2 + P2k']_ 



l£ 



2(?-2:b)(C-2:) [^2- 



kl 



j2 {^B-y){y-i) 
xb((,-xb) 



All integrals in k± we encounter during the calculation have the following form: 



7C-2; 



(27r)^-2 [^2 _ -g] [^2 _ JV^2 _ 

(27r)^-2 [fc2^ - Mf - is] [kl - Mi - le] 



a{M 



,r2\ 



a + PAP 

a + PMl 
Mf - Mi 
1 



a{M^) , 
a{Ml) - 



a + /3Af| 
Mf - Mi 



a{Ml) , 



(47r)- 



-M^ 



le] 



A-D 



so after k± integration wc get the following result: 



Att \47r/i2y {q + xp)^ + ie V 2/x — ^ 



2 a;s — X 



XB 



2 r^'^B 



dy 



y-x 



Xb — X 



i+f 



XB -y 
Xb — X 



Xb 



— le 



dy 



y-C 

Xb 



XB -y 
Xb 



M2 



-A 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 
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The last integration is performed making use of the beta function definition. The diagram with a left vertex correction 
is given by symmetry ^ — — ^, and the crossed diagrams hy xb ^ —xb- When we include all four vertex corrections: 



Tr [M'^'' f\ = ( Tr 



+ ^ -0 + {xB ^ -xb) 



(28) 



we get a result with the fohowiiig structure 
111 

= 5t 



^e^asCp ( 1 



47r \ 47r/i2 



Q 



.Xb — X 
Xb 



— le 



12 x-xe 



le 



le 



Q.X_B^ 

Xb 

{xb ^ -Xb) \ , 



fi (xb , X, e, e) + /i (xb , x,e,e) 

1 



^ Xb 



f2ixB,^,x,e,s) + 



^ SB 



f2i-XB,-^,x,e,e) 



(29) 



where /i and /2 are some complicated functions of xs, ^, x, e, 



C. Integrals with four propagators. 

All of the integrals with four propagators may be reduced to the three propagator case, although some of them 
require some care. We will start this section with the calculation of the symmetric part of the box diagram shown in 



Figure 4: Box diagram correction to q'y — >■ 57 scaterring amplitude 

the Fig. m 

^ 7 ' ' [(A: - $p)2 + le] [{k + ^p)2 + le] [{k + g)2 + le] [{k - xp)^ + le] ^ ' 

In this case we have four denominators in the integrated function, but one can easily check that: 

gr^, Tr C #)7''(^+ ^)ri^ + ^ lf)lp f{\ = A ■ k' + B ■ 2kp , (31) 

so, using following relations: 

- i(fc + ep)' + (c^-e), 

2k -p = l(fc + ^p)2 + (e^-0, (32) 
one can easily decompose the four denominator integral into two integrals with three denominators: 

.T,.Tr[A^-^ - iik - xp)^ + lik^-''S+ ^e] [(k + ,)2 + .s] + ^ ' 
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Figure 5: First group of diagrams describing 'yg "fg scattering. 



which wc calculate in the same way as the vertex corrections diagrams. Crossed diagram is given by the xb to ~xb 
replacement. 

Let us now turn to the gluon coefficient functions. The symmetric part of the first diagram describing 75 — ?> 
scattering, shown on the Fig. [5l is given by: 

CP,, _ . 2 2^ f,„^, 9T^p9T,. ■ Tr [7"(^ - X ^h^fi - ^ 4)1" + C #)] 

5t 9t A1(i)M.a0 -leg If j (dk) _ ^^^^ _^ _ ^^^^ ^ ^ ^ ^ ^^^^ ^ , [66) 

with Tp ~ i. The structure of the numerator is similar to the one given by equation (j3ip . So we can use the same 
decomposition as in the case of the quark box diagram. Diagrams (2), (3) and (4) from Fig. [5] arc connected to 
diagram (1) by simple symmetries. To get diagram (2) one has to change xb ^ —xb, diagram (3) x ~x, and to 
get diagram (4) one has to do both changes. 

Diagrams shown on Fig. [6l have different denominator structure, so we will describe the way of dealing with them 
more precisly. Momenta fiowing in the diagram (5) may be chosen as: 

A = k + q-{x-Op, 
B = k + q, 




Figure 6: Second group of diagrams describing 73 — >■ 717 scattering. 
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C = k + ^p, 

D = k-xp, (34) 
and Ag, Bq,Cq, Dq from diagram (6) are equal to A, B, C, D with q o — g. Both diagrams give the same result: 



, gTM..gTa;3Tr [7^ ^l" <PT P\ 

[A^ + ie\ [B^ + le] [C^ + ie] [D^ + le] 

As previously we notice that the numerator may be written &s A k^ + B 2k ■ p . To reduce our integral to the three 
denominator case, we use other relations: 

fc' = -^Ak + ipf + ^Ak-xp)\ 



2k- p = ——{k + ^pY —{k^xp) 

x + E, x + i 



so we end up with: 



1 Ax + B 1 r^,, A^-B 

^5) = TT-TT / W r42^,virR2^,virn2^,.i + TT-TY / (^^) 



= h+h. (36) 

One could worry if the above decomposition is well defined for x = — ^, but it is easy to check, that the expression 
(p6)) is regular in that limit. 

As previously we start with integration over k^ . We find four poles: 



(k - xpY + ie = 
{k + qY + ie = Q => k^ 
{k + q- {x- Op)^ + ie = ^ k^ 



le 



2(y-a2) 
kl -le + Q^f^ 



(fe ~ ^pf + ze = kl = , (37) 

2{y - 04) 

where ai = x, 02 = a^s, as = xb + x — ^ and 04 = — ^ are values of y for which poles imaginary parts change sign. 
Again appropriate choice of the integration contours allows us to write: 



^2 = ^ J Res,- /2 - ^, j^^ - j ^--^^ Res,- H • (38) 

The only difference with equations ([23]) is that we now have additional mass term in the denominator: 
Res,-/i 
Rcs,-/i 
Rcs,-/2 

r> J- ^ y "-4 "4 1 M4"'j_ /on^ 



1 






y~ai 






ai + /3ifc2 




a; + 


? 


2(ai 


- a2)(ai 


- as) 


[kl 


- Ml^ - ze] [kl - 


Ml, ~ ze] 


1 






y - ^2 






a2+p2kl 




x + 


? 


2(a2 


— ai)(a2 


- as) 


[kl 


- Mf, - ze] [kl - 


ie] 


1 






y-a2 






a3 + Pskl 




a; + 




2(a2 


- 04) (02 


- as) 


[kl 


- - ze] [kl - 


ie] 


1 






y - a4 






a4 + P^^kl 




x + 




2 (04 


- 02) (04 


- as) 


[kl 


- Ml - ze] [kl - 


A/I3 - ze] 
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where ai,f3i and Mij arc now defined by: 



ai,2 + Pi,2kj_ 
M,,, 



Ax + B 



q2 {y-at){aj -y) 



XsicLi - ttj) 



(40) 



Making use of Eq. ([26|) we arrive to 

1 



'(5) 



+ 



(47r)2 \4:Tr fi'^ 

a a 

ai (ai - a2)[ai - as) 
dy 



21 x + S, 



dy 



a4 
as 



(a2 - ai)(a2 - 03) 

y - Q4 

(04 - a2)(a4 - 03) 
y-a2 

ia2 (02 - a4)(a2 - 03) 
One can check that the relations: 

y-ai 



dy 



a2 + p2Ml^ 

Mh - Ml 
as + 



{-Ml, 
{-Ml, 



42 



Ml^ - Ml, 



a4 + /34MI3 
MI2 - Mi, 



-M^ 
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43 



42 



Ml 



42 



r 



(41) 



1 



(ai - 02) (ai - 03) A/f^ - Ml, 
y -a2 1 



(52 (ag - a2){y - ai) ' 
xb 1 



(02 - ai)(a2 - 03) MI2 (02 - a3){y - ai) ' 

ai + l3iMl = a2+P2Ml, 

allow to rearrange the integration limits, so we can express our integrals in the form allowing us to perform intej 
over y again using the beta function definition: 



(42) 

(43) 
^ration 



(5) 



i ( 1 



(47r)2 \^_'Kp?■ 
dy 

ai 

dy 

ai 
0.1 

dy 

as 

dy 



r 



y 


- ai 




1 


y 


- ai 




1 


y 


— 04 




1 


y 


— 04 



1 Xb 
27 a;2 _ ^2 Q2 

iai+^iMl){- Mi- 



le) 



le] 



{a,+P,Mt,){-Mt, 
(a4 + PaMD {-Ml - ie) 
iai + f3iMl){-M^ 



M3 



(44) 



In the next section we will write explicitly the final results of all of the above calculations. 



IV. RESULTS 



We see that hard scattering amplitudes for general kinematics have the following structure: 



e'^agCp 1 
47r (47r^2 



Q 



2X B - X 
Xb 



f{x,^,XB,e,£) 
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— le 



Xb 



(45) 



which in the e — > Umit for 57 — >■ (77 amphtude gives: 



Tr [M^"" f\ = 5^- 



47r 



12 



X — Xb + 

( i&{xxB - e) 



8(a: - Xb) 



{x-XB+ie^){x^ -e) 2)2-^2 



log ( sgn(Q^ 



.Xb - X 

) ■• 

Xb 

^ ix(xB - 

{x-XB+ie^){x-0 {x + XB -ie^){x + ^ ^(a^^ - C^) ^ 
18 



log ( .gn(g2)^^ 



X-XB+ 



-6 



a;2 +^2 _ 2xxB 



ix-XB+ie^){x^ ~e) 

' MxxB - e) 



log ( sgn{Q^)— — ^ - le 



2(a; - Xb) 



ix-XB+ie^)ix^ -e) x^~e ^ 

H^-xb) 6{^~xb) 



log' [sgn{Q')^^^^^e 



{x + XB ~ie^){x + i) (x - xb +is^){x - C) ^ 
2{^~xb) 2{^~xb) 



Xb 



log(s.9n(Q2).^^ ^ 



XB 



2x{xB - i) 



\{x-XB+ie^){x-0 [x + XB-ie^){x + i{x^-e) 
+ {xb < — > -Xb) , 

and for ^ gj- 



log^ sgn{Q') 



XB 

Xb 



(46) 



9t 9t ■'^tii^aP 



47r \47r/i2 

32{x^ -2xBX + 2xl -^'^) 



x^-e 

32ixB-Oix^ -2xB^-e) 

ax^ - e) 

-16 log ( sgn{Q^)^^^ - is 



log sgn{Q ) 



2,Xb~X 



l£ 



Xb 



log ( ,g„(Q2)^^ 
XB 



Xb 



(x^ - 2xbx + 2x% - ^2 
x2 -^2 

XB 



log' sgniQ') 



2 Xb — X 



— le 



Xb 



16 1 - 



log sgn{Q' 



Xb 



IS 



%(xB-i){x' -2xBi-e 



+ {xb 



ax^-e) 

-Xb) ■ 



log" sgn{Q' 



^ xb 
Xb 



(47) 



From the above result one can easily read the coefficient functions defined by Eq. ([7]), necessary for a calculation of 
the DDVCS amplitude, by means of Eq. (fT4|) . However, most experimental data is, or will be available for either 
DVCS or TCS, so we will elaborate more on those cases. 

Below wc present the resulting coefficient functions for the limiting cases of DVCS {Q^ > and xb ~ and TCS 

(-Q2 = Q'2 > and XB = -0- 
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A. DVCS limit 

We present the results explicitly showing ie terms that uniquely determine all imaginary parts. Quark coefficient 
functions read: 



1 



1 



47r 



a; — ^ + ie 



+ 



e^asCp ( 1 



4tt 



+ 



X — + ie 
1 



- £, — ie 

- 9 + 31og(l - ^ - ze) - 6^1og(l - ^ - ie) + 6-^log2 
C ^ x+^ 

+ log2(l_|_,£)_log2 2 

9 + 3 log(l + ^~ie) + 6^— log(l + I - le) - 6^- log 2 
+ log^(l + |-i£)-log^2 
6 + 41og(l - I - ie) -41og2 



a:; + ^ - i 
Gluon coefficient functions read: 



6 + 41og(l + - - ie) -41og2 



r<9 

'-'coll 



8x 



47r 



(x + ^ - ze)(a; - ^ + ie) 



^log (^1 - - -zej + ^log (^1 + ^ -.e 



2^^ log 2 



CI 



Att {x + ^ — ie){x — ^ + ie) 

X 



a; + 4 



a: 
? 

X 



le 



le 



(48) 



(49) 



Although the result for Cf contains dangerously looking denominators inside the square parenthesis, it is easy to 
check that the expression inside those parenthesis is regular in the limits x — > 

Above results are in agreement with earlier results [if 01 which were obtained in an unphysical region of parameter 
space, and then analytically continued to obtain DVCS case. We see that the simple prescripton that all imaginary 
parts can be obtained by substracting a small imaginary part from ^, i.e. ^ — >■ ^ — ie, is confirmed by our calculations. 



B. TCS limit 



Quark coefficient functions read: 

1 1 



^1 - 



''ya; — ^ — ie a; + ^ + ie 
1 



a: — ^ — ie 



-9 + 3 log(-l + ^ - ze) - 6^ log(-l + ^ - ie) + 6^ log(-2 - ze) 
i x + i £, .a: + 4 

+ log2(-l + I - ie) - log2(-2 - ze) 



13 



■ le 



9 + 31og(-l 



log'(-l 



X 

? 



ie) + 6- 



'-'coll 



An 



£, — ie 



41og(-l 
f 41og(-l 



X 



X 

1 



— ie) 
ie) - 41og(-2 

- ie) - 41og(-2 



log(-l 

— ie) 
■le) 

- ie) 



^ - ie) - 6—^ log(-2 - ie) 



(50) 



Gluon coefficient functions read: 

(Eg e?) as7> 



r<9 

^coll 



8x 



47r 



(a; + ^ + ie)(x - ^ - ie) 



x + S, 



log 



te 



log 



2x 



- 2 



- 2- 



3e 



a; 
. X - 



X - 

^2 



log 
log 



+ ie){x 

X 

1 



— ~ le 



X 



— 1 — — — ie 



log(-2 - ie) - 2- 



— ie) 



X 

? 



log^ -1 + - - ie 



■ log(- 



loa 



— le 



log^(-2-^e) 



(51) 



As in the DVCS case terms inside the square parenthesis of Cf are regular in the limits x 

There arc some important difFcrcnccs between the Eqs. ([50l [51]) describing the TCS case and Eqs. 
describing DVCS. First we notice that we have to add small imaginary part to ^, i.e. ^ — >■ ^ + ie, rather then substract 
as in the DVCS case. The second difference is the minus sign under the logarithms, which produces additional terms. 
Particularly log^(— 2 — ie) present in the TCS result may produce correction much bigger then the log^(2) in the 
DVCS case. Another important difference between the DVCS and TCS amplitudes concerns their imaginary parts, 
which in the DVCS case is present only in the DGLAP region, while in the TCS case, it is present in both DGLAP 
and ERBL regions. 



V. CONCLUSIONS 



It is well known that at the Born level TCS and DVCS hard scattering amplitudes are related: 



^O(DVCS) ~ ^O(TCS) 



The same relation holds for the coUinear terms: 



'-'colliDVCS) ~ ^coll{TCS) ' 



(52) 



(53) 



as they are equal to the convolution of the same evolution kernel with Born level amplitudes. Indeed this equality is 
crucial for factorization to hold. But in the NLO this relation no longer holds. For the quark part, we have : 



^1{TCS) ^l(DVCS) 



1 



X — ^ + ie 
1 



3-21og2 + 21og|l--| ) (i7r)+7r2(l + 6l(x-e) 



-^ + 0) 



3-21og2 + 21og|l + || ) (i7r)+7r2(l + 6l(-a;-e)-6'(a; + C)) 



a; + ^ — ie 

To discuss this difference and present the magnitude of corrections we define the following ratio: 



(54) 




Figure 7: Real (solid line) and imaginary (dashed line) part of the ratio R'' of the NLO quark coefficient function to the Born 
term in Timelike Compton Scattering (up) and Deeply Virtual Compton Scattering (down) as a function of x in the ERBL 
(left) and DGLAP (right) region for ^ = 0.3, for fi% = \Q'^\. 



coll 



R' = (55) 



^0 

of the NLO quark correction to the coefficient function, to the Born level. On Fig. [7] we show for = \Q'^\ the 
real and imaginary parts of the ratio R'^ in timelike and spacelike Compton Scattering as a function of x in the 
ERBL (left) and DGLAP (right) region for ^ = 0.3. We fix = 0.25 and restrict the plots to the positive x region, 
as the coefficient functions are antisymmetric in that variable. We see that in the TCS case, the imaginary part 
of the amplitude is present in both the ERBL and DGLAP regions, contrarily to the DVCS case, where it exists 
only in the DGLAP region. The magnitude of these NLO coefficient functions is not negligible. We see that the 
importance of these NLO coefficient functions is magnified when we consider the difference of the coefficient functions 
^i(TCS)* ~ ^i(DVCS)- '^'^^ conclusion is that extracting the universal GPDs from both TCS and DVCS reactions 
requires much care. 

As is well known in inclusive reactions, one may choose a renormalization scheme (named DIS scheme [l5| ) defined 
by the fact that NLO corrections to some observables vanish. This of course docs not preclude the importance of 
NNLO corrections. In the exclusive case, we thus may propose that NLO corrections vanish in the DVCS amplitude. 
This DVCS factorization scheme then transfers all NLO corrections calculated here to the TCS coefficient functions, 
which become very sizeable. We illustrate this fact by showing on Fig. [5] the ratio R^-s the difference of NLO 
quark coefficient functions to the LO coefficient function 

A final word is needed with respect to the presence of the tt^ terms in the difference of the NLO coefficient 
functions. Quite a rich literature [isl . [l6| exists on the importance of such factors in inclusive coefficient functions 
and their relation to soft gluon exchange. One may verify that in the exclusive case that we study here, a soft gluon 
approximation gives some of the tt^ terms that one may read from Eq. (|54)) . One can suppose that these corrections 
exponentiate when all order corrections are summed up. A particular feature is worth to be pointed out : these tt^ 
terms only exist in the DGLAP regions. We confess that we do not understand why this is the case. 
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Figure 8: Real (solid line) and imaginary (dashed line) part of the ratio R^_g of difference of NLO quark coefficient functions 
to the LO coefficient functions in the TCS and DVCS as a function of x in the DGLAP region for ^ — 0.3. 
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Figure 9: Ratio of the real (solid line) and imaginary (dashed line) part of the NLO gluon coefficient function in TCS to the 
same quantity in DVCS as a function of x in the DGLAP region for ^ = 0.05 for fj,% — \Q'^\. 



Let us now briefly comment on the gluon coefficient functions, 
parts arc complex conjugated to eacli other: 



As in the case of quark corrections, the collinear 



^coU(DVCS) 



^coU{TCS) 



(57) 



Moreover, the real parts of the gluon contribution are equal for DVCS and TCS in the ERBL region. The differences 
between TCS and DVCS emerges in the ERBL region through the imaginary part of the coefficient function which is 
non zero only for the TCS case and is of the order of the real part. In the DGLAP region, the difference reads : 



^1{TCS) 



^1{DVCS) x>i 



2x 



47r 



x + £, 



log|l-|| 



(58) 



showing a sizeable difference of the contributions to both the real and imaginary parts of the amplitude. In Fig. [9] 
we illustrate the ratio of the NLO gluon correction to the hard scattering amplitude in TCS to the same quantity in 
the DVCS in the DGLAP region for ^ = 0.05 for fj.^ = \Q'^\. 

The discussion of NLO corrections to a hard scattering amplitude necessarily brings the question of the factorization 
scale dependence. On the Fig. [TO] we show the real and imaginary parts of the ratio i?"? of NLO quark correction to 
hard scattering amplitudes to Born level coefficient function of the timelike Compton scattering as a function of x in 
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the DGLAP region for ^ = 0.05. The figures are plotted for various values of -t^, and present a strong factorization 
scale dependence. 

On the Fig. [TT]we show the ratios of the real (left) and imaginary (right) parts of NLO gluon coefScient function 

2 

for IQ^I = ^ (solid line) and |Q^| = 2yu|, (dashed line) to the same quantities with \Q'^\ = iJ?p. Those quantities are 
calculated for the timelike Compton scattering and plotted as a function of x in the DGLAP region for ^ = 0.05. Also 
in this case we notice a strong factorization scale dependence. 

Much phenomenological studies need now to be performed, by convoluting the coefficient functions to realistic 
GPDs and calculating the relevant observables. We plan to progress on these points in the near future. 




Figure 10: Factorization scale dependence of the real (left) and imaginary (right) parts of ratio i?' of NLO quark correction 
to hard scattering amplitudes to Born level coefficient function of the Timelike Compton Scattering as a function of x in the 
DGLAP region for ^ = 0.05. The ratios are plotted for the values of equal 0.5 (dashed), 1 (solid) and 2 (dash-dotted line). 




Figure 11: Ratios of the real (left) and imaginary (right) parts of NLO gluon coefficient function for IQ'^I = l/2/^|. (solid line) 
and IQ'^I = 2/i|. (dashed line) to the same quantities with |Q^I = A- Those quantities are calculated for the timelike Compton 
scattering and plotted as a function of x in the DGLAP region for ^ = 0.05. 
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